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We consider the coherent stimulated Raman process developing in an optically dense disordered 
atomic medium, which can also incoherently scatter the light over all outward directions. The Raman 
process is discussed in context of a quantum memory scheme and we point out the difference in its 
physical nature from a similar but not identical protocol based on effect of electromagnetically 
induced transparency (EIT). We show that the Raman and EIT memory schemes do not compete 
but complement one another and each of them has certain advantages in the area of its applicability. 
We include in our consideration analysis of the transient processes associated with switching off/on 
the control pulse and follow how they modify the probe pulse dynamics on the retrieval stage of the 
memory protocol. 
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I. INTRODUCTION 

Quantum memories for light pose an extremely impor- 
tant problem for various protocols of quantum computing 
and secure communication [l|, |2|, |3| , which might be solv- 
able with currently existing technology. At present there 
have been several experimental demonstrations of storage 
and retrieval of non-classical states of light from macro- 
scopic atomic systems 0, i, fl 0, H, !, El El El El ■ In 
many theoretical discussions the problem was tested in 
context of achievable fidelity for storage and retrieval of 
the original state of light, normally assumed either co- 
herent, squeezed or single photon 0, H H, EE EE El, E3 • 

In the recently reported series of papers [3] thorough 
analysis was applied to determine the efficiency of stor- 
age and retrieval of a single mode quantum state of light 
via its lambda-type conversion into an atomic spin coher- 
ence. Such a conversion is normally associated with ei- 
ther electromagnetically induced transparency (EIT) two 
photon resonance, or with Raman process, or with pho- 
ton echo effect. Optimization procedure and practical 
recommendations are often based on the unique proper- 
ties of the lambda transition and presumes the quantum 
information encoded in collective variables of the light. 
For the weak coherent light (with small degeneracy pa- 
rameter) this scheme describes, for example, the storage 
of a single photon quantum bit and relates the efficiency 
of the total pulse retrieval to the fidelity of the single pho- 
ton quantum state retrieval. Such analysis is relevant for 
one mode memory channel and it leaves open important 
issues concerning storage and retrieval of a multimode 
quantum state. The quantum memory for the multimode 
configuration in the Raman process has been the subject 
of recent reports [3, E3 , which indicated the difficulties 
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and necessity of additional studies. 

In the present paper we shall discuss the stimulated 
Raman process in a scattering medium with a A-type 
excitation, which is initiated by a control pulse of lim- 
ited duration. We shall split this mechanism from the 
similar but not identical EIT two-photon resonance and 
discuss the difference between these processes. The Ra- 
man scattering in the Heisenberg formalism, introduced 
in Ref.[3|, and later applied to the quantum memory 
problem in Refs. [3, EE E3 , is usually discussed for far 
frequency offsets of the probe and coupling modes from 
an atomic resonance and it is considered in approxima- 
tion based on adiabatic elimination of the upper states. 
This approximation makes possible to build the relevant 
effective Hamiltonian and include the losses into the dis- 
sipation terms and the Langevin forces. In the present 
paper we show how the Raman process can be described 
for arbitrary frequency detunings and without such a sim- 
plifying approximation. We show how the Raman process 
can be discriminated from the EIT two-photon resonance 
for any frequency detunings via considering the spectral 
behavior of the sample susceptibility dressed by the cou- 
pling field. 

Our analysis aims to clarify the physical requirements, 
for dependable manipulation with the light pulse and 
atomic coherence at each step of the complete memory 
protocol. The important physical processes like atomic 
motion and transient effects associated with switching 
off/on the control pulse will be described. In the last part 
of the paper we briefly comment how our approach can 
be generalized to a multilevel configuration while pay- 
ing attention to the hyperfine structure of alkali atoms, 
which was earlier considered in a lossless configuration 
m Ref.[l7l. As an important consequence of our anal- 
ysis we shall point out that the EIT resonance and the 
Raman process are not competitive protocols but actu- 
ally complement one another. The protocol based on the 
EIT effect works better for pulses longer than the natural 
atomic lifetime 7 , but the Raman protocol is in princi- 
ple applicable for pulses with duration around 7 -1 . Thus 
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the Raman scheme can potentially work for pulses around 
20 ns and hence extend the range of pulses, which could 
be stored in the atomic spin subsystem. That would lead 
us to wider demonstrations of atomic memory for un- 
known states of light arriving from radiation sources of 
different physical nature. 

The paper is organized as follows. In section |TT] we 
describe our approach to the stimulated Raman process, 
which is based on the Green's function formalism. We 
indicate some technical advantages in applying diagram 
analysis particularly to the atomic memory problem. In 
section IIIII we make a diagram derivation of the master 
equation describing the dynamics of the probe pulse and 
atomic coherence on the write-in and retrieval steps of 
the memory protocol. Our general discussion is further 
illustrated by numerical simulations for the Raman pro- 
cess given in section IIVI Some details concerning the 
Green's function formalism are presented in appendix [A"l 

II. BASIC ASSUMPTIONS AND 
CALCULATION APPROACH 

Transport of a light pulse in a medium under condi- 
tions, where the coherent forward propagation is damped 
by losses due to incoherent scattering over all outward di- 
rections, is a quite delicate problem to describe theoret- 
ically at the quantum level. The problem becomes even 
more subtle for practically important situation when the 
energy of the transmitted pulse is comparable with the 
scattered light energy. If the losses are significant and 
the light and atomic subsystems are initially prepared in 
arbitrary quantum states it would be difficult to follow 
the dynamics of the entire atoms-field subsystem com- 
pletely in a three dimensional configuration, see [2C| . An 
essential simplification can be achieved if the light and 
atomic subsystems are initially in a coherent state and 
the interaction process does not modify this type of state 
[2ll ] . Such a dynamical process, which reproduces most of 
the principle features of the quantum memory protocol, 
we shall discuss in this section. 



A. Description of the light subsystem 

For a coherent state any quantum characteristics of 
light can be defined in terms of the complex amplitude 
of the electric field given by the expectation value of 
its Heisenberg operator considered at the certain spatial 
point r, time t and averaged over the asymptotic state of 
the system in the infinite past 

4+>(r,t) = (4 +) (r,i)> = ($TE${r,t)$) (2.1) 

Here £^ + ' (r, i) and Squ ( r > are tne positive frequency 
parts of the /i-th polarization component of the electric 
field operators in the Heisenberg and interaction repre- 
sentations respectively, which are related by the following 



unitary transform 

4+>(r,i) = St(i,-oo)4?M)S(i,-oo) 

S(t,f) = Texp(-^J V {t")dt") (2.2) 

where Vo(t) is the interaction Hamiltonian in interac- 
tion representation, see Eq. (|2.6p below, T is the time- 
ordering operator, and the S-matrix in expression (|2.ip 
is given by principle limit of the evolutionary operator 
S = S(oo, —00). 

If the light pulse crosses the medium in a single pass, its 
forward propagation obeys the mesoscopically averaged 
Maxwell equation, which can be constituted and solved 
under rather general assumptions. The field emerging 
from the sample via incoherent scattering can be visual- 
ized as secondary waves created by a disordered system 
of atomic sources. For the weak field, scattered by the 
sample, its intensity, angular and spectral distribution 
can be found by analysis of the multiple scattering pro- 
cess. To entirely simulate this process the Monte-Carlo 
scheme can be applied. 

Let us make these general statements more concrete 
in context of the A-configured transition, shown in fig- 
ures Q] and [21 The depicted scheme is routinely discussed 
as theoretical background for description of the stimu- 
lated Raman process or EIT two-photon resonance. A 
strong coupling field with right-hand circular polarization 
is quasi-resonantly applied between a Zeeman hyperfine 
sublevel of the ground state \m') and a hyperfine com- 
ponent of the excited level near the upper Zeeman state 
|n). A weak probe beam is also quasi-resonant with the 
\m) — > |n) Zeeman hyperfine transition. All nonlinear op- 
tical effects associated with the probe light are assumed 
to be negligible. In our discussion this field will be taken 
into account only in the first non-vanishing order. The 
atomic medium consists of the atoms initially (in the in- 
finite past) populating state |m), such that the presence 
of the weak probe initiates only weak transfer of atoms 
from state \m) to \m'). The sample is optically dense 
for the probe field and fully transparent for the coupling 
field. 

Both the coupling and the probe fields are in a coherent 
state and completely characterized by the complex ampli- 
tude (|2.ip projected onto the relevant polarization and/or 
spectral states that makes these modes distinguishable. 
The properties of the coupling field are not modified at 
all because the sample is transparent and the number of 
photons coherently scattered into the coupling mode is 
negligible in comparison with the quantum uncertainty of 
the number of photons contributing to this mode. Thus 
the evaluation of the expectation value (|2.ip based on 
a relevant diagram expansion need to be performed only 
for the probe field. This calculation will automatically in- 
clude all the dressing effects associated with the coupling 
field and describe in one formalism two different physi- 
cal phenomena: stimulated coherent Raman process and 
EIT two-photon resonance. The difference between these 
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FIG. 1: (Color online) Geometry of the quantum mem- 
ory scheme based on stimulated Raman process. The left- 
handed polarized probe pulse with carrier frequency Gj is co- 
herently scattered into the right-handed polarized strong cou- 
pling mode ui. The action of the coupling field is limited in 
time how it is indicated by the shadowed area in the temporal 
dependence of the process for input/output probe light inten- 
sity I(t). Thus the part of the probe light will not emerge 
the sample and will be transformed into the long-lived spin 
coherence in atomic subsystem. 



phenomena will be further clarified by considering differ- 
ent physical limits and spectral relations attaining the 
delay effect and minimizing the incoherent losses. 

The global advantage of the coherent state in compar- 
ison with other more sophisticated quantum states, such 
as squeezed, entangled etc., is a relatively simple diagram 
analysis for the expectation value of the field amplitude 
(|2.1j) . Its expansion in the series of perturbation theory 
under application of Wick's theorem transforms any out- 
coming product of the field operators to iV-ordered form, 
and generate only the simplest diagram objects such as 
the non-perturbed coherent amplitudes of incoming field 
and the vacuum-type Green's functions. This makes it 
technically possible to further describe the probe field 
dynamics in terms of mesoscopically averaged Maxwell 
equations and multiple scattering light diffusion. 



B. The atomic subsystem and interaction process 

We consider the macroscopic atomic system in the sec- 
ond quantized formalism. In this subsection we shall 
show that for an atomic ensemble underwent coherent- 
type interaction process shown in figures [T] and [2] the 
quantum dynamics can be properly described by the fol- 
lowing four Green's functions of the Perel-Keldysh dia- 
gram technique [H [H, Hi] 



,, CT * m ,(r',0*lM)) (2.3) 



These Green's functions with cr, a' = =F give the expecta- 
tion values of the different products of creation and anni- 
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FIG. 2: (Color online) Energy diagram illustrating in exam- 
ple of a A-configured transition how the dressing effects mod- 
ify the probe pulse propagation under the stimulated Raman 
process, shown in figure [T] The probe pulse with a carrier 
frequency u is scattered into the strong mode of frequency u> 
and creates coherence between the lower Zeeman states \m) 
and \m'). In the next step the Zeeman coherence initiates the 
delayed retrieval of the probe pulse. If the coupling mode is 
"cw" with infinite pulse duration then the propagation of the 
probe pulse under the Raman process is damped by the spon- 
taneous scattering on the components of the Autler-Townes 
doublet shown in red-dashed. 



hilation ^-operators ordered by chronological operators 
T a i a . The permutation operators realize the following 

action: T = T arranges the ^-operators as ordered in 

time, its counterpart T++ = T arranges them as anti- 

ordered in time, operators T_| and T |_ are respectively 

identity and transposition operators. 

The ^-operators responsible for anihilation or creation 
of an atom at spatial point r and in internal state \m) 
can be expanded in the basis set of plane waves as follows 



* m (r,t) = £ J=e^ r 6 pm (i) 
p 

*l(r,t) = ^J=e-ib r &t ro(t) (24) 



where V is the quantization volume, and b prn (t) and 
& pm (i) are respectively the operators of annihilation or 
creation of an atom with momentum p in internal state 
\m) in the Heisenberg picture. 

If the atomic motion has semiclassical description then 
the Green's function G^,^(r', t'\ r, t) before interaction 
can be conveniently expressed by the single particle den- 



4 



sity matrix in the Wigner representation 



[2] is driven by the following interaction Hamiltonian 



^LTLV,i';r,t)=± 



d 3 p 
(2irky 



exp 



-p(r' - r) 



x exp 



xp 



mm \ ft 



if-( t > - t )- i E m'+E m , _ 

h2m y ' h 2 V ' 
r' + r t' + 1 



(2.5) 



where the external sign + or — corresponds to bosonic 
or fermionic statistics of atoms. In this integral rela- 
tion the exponential factors are responsible for free dy- 
namics of the atomic state and show fast oscillation in 
space and time. The Wigner density matrix p m 'm(p, r, t) 
is also oscillating in time with the frequency u) m > m — 
(E m / — E m )/Ti such that the integrand itself has free pre- 
cession as function of times t and t' with internal energies 
E m and E m i respectively. The validity of expression (|2.5|) 
is constrained by the requirement that the atomic de- 
Broglie wavelength A a ~ %/ v/mT, where "m" is atomic 
mass and T is ensemble temperature, should be much 
less than the scale of spatial inhomogeneity. Then the ap- 
proximation of the true Green's function G^T/i, (r', t'; r, t) 
by expression (|2.5p is applicable if variations of the in- 



crements r' — r and t' — t are small in comparison with 
the spatial and time scales where the Wigner function 
is significantly changed. For free dynamics and without 
spin relaxation process the density matrix of the ground 
state can only be modified because of atomic motion in 
free space or in a trap. 

It is important to recognize that without interaction 
with the light subsystem the Wigner density matrix 
Pm'm(P) r : t) gives a complete description of the macro- 
scopic atomic subsystem. Indeed, as one can verify 
each of the remaining Green's functions (|2.3[) can be ex- 
pressed by integral relations similar to expression (|2.5p . 
see [23|, [2J] for details. The additional terms, which dis- 
tinguish the four Green's functions, given by Eq. (|2.3[) . 
are expressed by the commutator (for fermion case by 
the anti-commutator) of the free Heisenberg ^-operators. 
These terms contribute to G^ ' with additional multi- 
plication by 9(t' — t) to G( _ ~) (retarded type) and by 
9{t — t') to G' ++ ) (advanced type). In turn each commu- 
tator (or anti-commutator) can be expressed by a simi- 
lar expansions in the momentum representation and de- 
scribed by similar exponential functions of r' — r and t'—t 
as in Eq. (|2.5p if the atomic motion is free. For a coherent 
initial state of the atomic system, when all the atoms oc- 
cupy only one Zeeman state, any correlation function of 
higher order can be expressed as a product of the lower 
order Green's functions defined by Eq. (|2.3[) . 

The coherent type interaction process shown in figure 



Mt) ="£/ d 3 r(dnn m E^\r,t)^ 0n (r,t)^ 0m (r,t) 



H.c. 



(2.6) 



where all the time dependent operators are written in 
the interaction representation as indicated by the zero 
subscript and (d tJ ') nm denotes the matrix element of the 
/U-th component of an atomic dipole. The second line 
constitutes the rotating wave approximation and convert 
to algebraic form the normal differential relation between 
the transverse vector potential operator Aqu (r, t) and 

the electric field operator E^\r,t). This presumes that 
both the frequencies of the coupling field u> and the probe 
field a) are indistinguishable from an average frequency 
too for the atomic optical transitions shown in figure [2] 

As one can verify via perturbation theory to lowest 
order the probe field, the interaction process shown in 
figure [2] will not break the Gaussian factorization for any 
higher order correlation functions. This means that the 
interaction process leaves the atomic system in the state 
described by a single particle density matrix. Any atomic 
correlation function is described by the retarded or ad- 
vanced type Green's functions dressed by interaction with 
the strong field. As we shall further show the main op- 
tical characteristics of the sample are expressed by such 
a dressed retarded Green's function of the excited state. 
This Green's function can be found in analytical form for 
a control pulse of arbitrary duration and in quite general 
assumptions, see appendix [A} The propagation of the 
probe pulse through the medium can be visualized as a 
polaritonic wave, both for the stimulated Raman and the 
EIT processes. 



III. DYNAMICS OF THE PROBE FIELD AND 
ATOMS FOR A-TYPE EXCITATION IN THE 
OPTICALLY THICK SAMPLE 



A. Write-in stage of the memory protocol 



The dynamics of the entire system consisting of the 
probe field and atoms can be described via diagram anal- 
ysis of the atomic and field Green's functions introduced 
in the previous section. After a certain regrouping of par- 
tial contributions generated by the expansion of pertur- 
bation theory for the evolutionary operator one obtains 
the following Dyson equation for the probe field ampli- 



5 



tude 




m 



(3.1) 

Here the double wavy single ended lines describe the 
probe field amplitude dressed by the coherent interac- 
tion shown in figure [2] The single double ended wavy 
line is the causal (retarded in the RWA approach) Green's 
function for light propagating freely in vacuum. The loop 
consisted of the atomic Green's functions and the vertices 
describe the polarization operator or susceptibility of the 
sample in response to the probe field. The equation (|3.ip 
graphically represents the macroscopic Maxwell equa- 
tion. More detailed analytical description of the atomic 
Green's functions is specified in appendix IA1 

The dynamics of the atomic coherence can be identified 
if the solution of equation (|3.ip is known. This requires 
evaluation of the following diagram for the Green's func- 
tion G ( ~,X}{- ■ ■) defined by Eq. (j23]) and created by the 
interaction process 



The Maxwell equation (|3.3p indicates that in a dilute 
sample the coherent part of the probe pulse propa- 
gates only forward and the profile of the outgoing pulse 
strongly depends on properties of the sample suscepti- 
bility x(r±,z;t,t'), which controls the efficiency of the 
memory protocol. The diffraction divergence was com- 
pletely ignored in the derivation of Eq. (|3.3|) and r j_ was 
treated only as a parameter. 

Since the coherent process shown in figure does not 
disturb the original equilibrium distribution of atoms, the 
ground state density matrix can be factorized as 



(p,r,i) = a m 'm(r±,z,t) / (p) 



(3.5) 



where /o(p) is the equilibrium Maxwell distribution of 
the atomic momenta. In the case of an isolated single 
A-type excitation the diagram expression (|3.2[) leads to 
the following kinetic equation for the off-diagonal ground 
state matrix elements of the atomic density matrix 



0_ 

at 



(r ± ,z,t) + (r ± ,z,t) =A£,, 



xe -i(u,- u )t+i{k-k)z dt'x(r ± ,z;t,t')e(r ± ,z;t') (3.6) 



where the scaling factor 



m 



" ' + (3.2) 

Here the internal doubly line represent the retarded type 
Green's function of an excited atomic state dressed by 
the interaction with the vacuum and coupling modes, see 
appendix [A] The outward directed dashed arrow repre- 
sents the field amplitude of the coupling mode. The thin 
lines are related to the non-perturbed dynamics of free 
atoms in the ground state and are given by semiclassical 
approximation (|2.5[) . 

The strategic idea of the write-in step of the mem- 
ory protocol is that equation (|3.1[) should generate probe 
pulse delayed so much as it practically does not emerge 
from the sample during the action of the control pulse. 
At the same time the probe light should not be incoher- 
ently scattered by the sample and after the interaction 
the information, originally encoded in the probe pulse, 
will be concentrated in the spin subsystem. The dynam- 
ics of the probe field in this process is driven by equation 
(|3.ip . which has the following analytical form 



Id d 1 . . u> f , . _ . ... i. 

cdt + dz e y r - s -' z '^ t ) =2m ~ / dt X{ r ±,z;t,t)e(r ± ,z;t ) 

— OC 

(3.3) 

where e(r±,z;t) is the slow- varying amplitude of the 
probe mode, which is given by a factorization of the probe 
field component defined by Eq. (|2.ip 



4rftM)=e(r±,z;i)e 



V* 



(3.7) 



is responsible for the difference in the interaction vertices 
contributing to the polarization and self-energy opera- 
tors of the diagrams (|3.ip and (|3.2p respectively. Here 
V nm ' = (de_|_) nm /£ is the transition matrix element for 
interaction with the strong coupling field, and \m'), \m) 
are the Zeeman states coupled by the interaction process. 
Vectors e+ and e_ are respectively unit polarization vec- 
tors for right-handed and left-handed polarizations. The 
strong mode is assumed to be described by the following 
positive frequency component 



42 h t(M)=£e- 



S(t - z/c) 



(3.8) 



i. e. to be nearly monochromatic during the action of 
the control pulse with the envelope profile S(t). 

Both equations (|3.3[) and (|3.6p are described by the 
same sample susceptibility. The hardest obstacle for fur- 
ther application of equation (|3 .3[) is that the susceptibil- 
ity is not known analytically in the general case, when 
atoms have a multilevel energy structure and the con- 
trol pulse have an arbitrary envelope profile. It can be 
described analytically only after a round of approxima- 
tions with respect to the dressing effects associated with 
the coupling field and the internal atomic interaction. In 
appendix [A] we calculate it for the case of rectangular 
profile, when the control pulse has duration T 



-iujt+ikz 



(3.4) 



S(j) = 0{t) - 9(t - T) 



(3.9) 
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where 0(t) is the step function. For this particular case 
the Maxwell equation (|3.3[) can be solved numerically. 
Then the distribution of the spin coherence in space and 
in time can be found via straightforward integration of 
equation (|3.6p . 

Together the equations (|3.3|) and (13. 6|) describe trans- 
port of a spin polariton mode through the atomic sample, 
which dynamics can be approximated by either EIT-type 
[IH or Raman-type As was mentioned above, the 

strategic line of the write-in step of the memory proto- 
col is to stop the control pulse at a time, when the probe 
mode has not emerged from the sample and the polariton 
is localized at least on the scale of the sample. Then after 
a round of transient processes the state of the light will 
be mapped onto spatially extended state of the atomic 
spins in the sample. Without any external disturbances 
this state will be preserved and will store the mapped 
quantum state of the probe light. 

In ideal situation this process does not add extra 
noise to the state. The crucial requirement to make 
the memory protocol feasible, as it was discussed in 
Refs.0, [H, Hil, is that the sample should have ex- 
tremely high optical thickness, such that noA 2 L ^> 1, 
where no is a density of atoms, A is a wavelength di- 
vided by 2ir for either probe or coupling mode, and L is 



the sample length. If this requirement is fulfilled, which 
is principally important for both the Raman and EIT 
memory schemes, the delay time for the transport of the 
probe pulse through the sample can be made longer than 
the input pulse duration and the spin polariton can be 
localized in the sample. 



B. Retrieval stage of the memory protocol 

If the coupling light is switched off and after a con- 
trollable delay switched on again, the probe pulse can be 
stored and retrieved. The delay time is mainly limited by 
relaxation processes for the spin coherence. In the case 
of cold atoms the spin relaxation can be minimized but 
atoms can leave the interaction area because of free mo- 
tion while the confining trap fields are switched off dur- 
ing the entire protocol time. However for delay up to few 
milliseconds the atoms can be considered as non-drifting 
and the retrieval stage can be initiated by recovering the 
polariton dynamics with the second coherent pulse. 

The probe field dynamics is now described by the fol- 
lowing Dyson equation 





(3.10) 



r 



The first term on the right hand side displays the source 
for the optical coherence initiated by the second control 
pulse. The black wavy line performs here electric field of 
the incoming probe pulse, i. e. the solution of equation 
(|3.ip . Since the probe pulse has arrived with the write- 
in control pulse the right vertex of this term relates to 
times associated with the write-in stage of the protocol. 
But the left vertex relates to times belonging the readout 
control pulse. 

The broken double line denotes the retarded Green's 
function of the excited atomic state dressed by both the 
write-in and readout control pulses and has the following 
diagram definition 



tex points associated with the coupling fields for write-in 
(W) or readout (R) control pulses. Just because the vac- 
uum dynamics, described by the thin solid line, preserves 
the atomic state for infinitely long time this diagram be- 
ing inserted into equation (|3.10p can be visualized as the 
spin coherence, which was originally created by the pro- 
cess (I3.2j) . and which now regenerates the weak probe 
field via a stimulated Raman process when the readout 
control pulse is applied. 



R 



W 



(3.11) 

The incoming double line is the excited state Green's 
function dressed by the write-in control pulse. The out- 
going double line is the excited state Green's function 
dressed by the readout control pulse. These two lines are 
converted into the ground state vacuum line at the ver- 



The second term in equation (|3.10[) is similar to the 
second term in equation (|3.1[) and is responsible for the 
Maxwell dynamics of the recovered probe field. The red 
dashed wavy line reproduces the recovered probe field for 
which this equation actually should be solved. In analyt- 
ical form the diagram equation (|3.10p can be written as 
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follows 



ld_ d_ 

c dt dz 



£out(i"_L,z;t) 

t 

2m— I dt'x // (r ± ,z;t,t')e(r± 7 z;t / ) 



-oo 

t 



+2m- dt'x(Tx,z;t,t')e ont (r x ,z;t') (3.12) 



where e(r±,z;t') is the solution of Eq. (|3.3p and 



£out(r±, z; t') is field amplitude of the recovered pulse, 
for which equation (|3.12|) should be solved. It is assumed 
that at the beginning of the readout control pulse the 
incoming probe pulse is over. Despite the lower limits of 
the integrals in the right hand side are in infinite past, in 
reality the integrands are non- vanishing only if t' belongs 
the writing control pulse in the first integral and the read- 
out pulse in the second. The susceptibilities X//(- • •) an d 
x(. . .) are given in appendix |A1 in example of an isolated 
A-type transition. 

The dynamics of atomic coherence is now described by 
the following diagram equation 





(3.13) 



r 



which in turn leads to the following kinetic equation for 
the density matrix (|3.5|) 



xe 



-i{uj— uj)t-\-i(k — k)z 



dt'x//{r ± ,z;t 7 t / )e{r J _,z;t') 



(3.14) 



Similar comments to those given for (|3.12p above apply to 
this equation, which is responsible for atomic dynamics at 
the retrieval stage. In application it is only necessary to 
solve equations (|3.3f) and (|3.12p . and then a m ' m (r±, z, t) 
can be found by straightforward integration of (|3.14[) . In 
practice it is even not so important to know the solution 
for remaining atomic coherence because all the retrieved 
information goes out with the recovered probe pulse and 
the equation (|3.14[) should lead to a trivial result in the 
infinite future cr m r m (r±, z,t — ► 00) — > 0. However we in- 
troduce this equation for completeness of the description 
of the memory protocol. 



IV. RESULTS AND DISCUSSION 
A. The EIT effect and Raman process 

Let us briefly discuss the difference between these two 
physical mechanisms when the propagation of the probe 
field is controlled by a stationary coupling field. Referring 
to the definitions of figure [2] one can define the following 



frequency detunings for the coupling mode A = u — u nm > 
and for the probe mode A = Co — u) nm from the atomic 
transition frequencies w nm i and u> nm respectively. Then 
one can select two important resonance relations between 
the carrier frequency of the probe mode u> and the energy 
spectrum of the "dressed" atom. The resonance point 
A = A defines the two photon resonance, which is asso- 
ciated with EIT. Indeed the propagation function of the 
atomic excited state (|A3|) as well as the sample suscep- 
tibility (|A13[) vanish at this point for a pure monochro- 
matic A-type interaction of the modes lo and Q. The 
sample becomes transparent and any spectrally narrow 
probe pulse distributed in the vicinity of this point will 
cross the sample with extremely low group velocity. The 
light delay, associated with the EIT effect, does not prin- 
cipally depend on frequency offset A and, in particular, 
the group velocity near the EIT point depends only on re- 
lation between the Rabi frequency of the coupling mode, 
atomic density and 7 and not on A, see Ref.[2||. For a 
sample with high optical depth, such that no A L 3> 1, 
and for any A the probe pulse characteristics can always 
be adjusted to realize its slow transport or storage in the 
sample. 

Another resonance situation takes place if the carrier 
frequency of the probe pulse is scanned near the fre- 
quency shifted quasi-energy component of the Autler- 
Townes doublet as shown in figure [21 This resonance 
point manifests itself as a narrow high amplitude spec- 
tral feature in the behavior of the dispersion and absorp- 
tion components of the sample susceptibility, see Ref. [26[ . 
In the Maxwell equation (|3.3[) the spectral overlap of 
the probe pulse with the dispersion part of the Autler- 
Townes resonance is responsible for lossless coherent for- 
ward propagation, which is normally associated with the 
stimulated Raman process. The overlap with the absorp- 
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tion part generates the incoherent losses caused by quan- 
tum radiation coupling of such an exciton-type atomic 
state with the continuum of vacuum modes. 

As can be verified in the limit when the Rabi frequency 
ft c = 2\V nm '\/h is much less than detuning A, the equa- 
tions (|3.3[) and (|3.6[) can be transformed to the system 
of equations introduced in Ref. and discussed later in 
the context of quantum memory in Refs. fbH. [l6l [l7j . This 
transformation is equivalent to the standard procedure 
of adiabatic elimination of the optical atomic coherence, 
which was a basic assumption of the above papers. But 
we can point out that such an approximation is only valid 
if in the integrand in Eq. (|3.3|) the carrier frequency of the 
probe mode is located near the Autler-Townes resonance 
peak, which always requires for large A some overlap of 
the pulse spectrum with this resonance. The importance 
of this overlap for the problem of optimization the probe 
pulse storage was earlier commented in Ref. fl8| and in 
particular we address the reader to the second paper 
of that series. It is also important that in the coher- 
ent stimulated Raman process the photons of the probe 
mode transform to the coupling mode in result of inter- 
action with all atoms of the ensemble. After a certain 
delay the photons of the probe mode are recovered via 
reversing interaction of atomic coherence with the cou- 
pling field. The Raman mechanism would exist even if 
the strong coupling field were exactly on resonance with 
the atomic transition, i. e. for A = when the problem 
becomes symmetric with respect to both components of 
the Autler-Townes doublet. We demonstrate this possi- 
bility by results of our numerical simulations given below. 

Applying the Raman process at any A the spontaneous 
losses can be always minimized by appropriate adjust- 
ment of the external parameters and by optimization of 
the pulse shapes, see Refs.(l6l. [Tsj . This normally as- 
sumes either optimization of the probe pulse parameters 
with respect to a given shape and the frequency detuning 
of the control pulse or alternatively the optimization of 
the control pulse with respect to a given probe pulse. In 
an idealized lossless situation the stimulated Raman pro- 
cess could be considered as evolving purely dynamically. 
But because of the losses the efficiency of the coherent 
forward scattering is not responsibly high. In contrast 
to the EIT transparency point, for the Raman situation 
the forward coherent scattering is more sensitive to the 
frequency offset A because the spectral behavior of the 
sample susceptibility near the Autler-Townes resonances 
is essentially modified with varying A. Still at a given 
detuning and high optical depth the parameters of the 
probe pulse can be adjusted for its optimal storage. 

Both the Raman and EIT mechanisms for delay of the 
probe pulse are generally applicable for any detuning of 
the coupling mode from the atomic resonance. As com- 
mented above, in a stationary situation, both are de- 
scribed by the spectral behavior of the same sample sus- 
ceptibility but considered in different spectral domains. 
For general relation among external parameters these 
spectral domains can be only conventionally discrimi- 



nated. However in many practical realizations the EIT 
zone and the Raman zone can be clear separated. It is 
important to recognize that for both the processes the 
crucial feature is the spectral behavior of the dispersion 
part of the sample susceptibility. It can be steep cither 
near the transparency point (EIT) or near the absorption 
point of the Autler-Townes resonance (Raman) . 

The delay effect attains by the spectral overlap of the 
pulse spectrum with the sample susceptibility near any 
of these points. For the EIT scheme it is preferable to 
keep the pulse spectrum closer to transparency window 
and the protocol would work better for longer pulses (in 
scale of natural lifetime) and for smaller Rabi frequencies 
of the coupling field. For Raman scheme the pulse spec- 
trum should be concentrated close but outward from the 
absorption peak. Thus the Raman protocol can be de- 
signed for large Rabi frequencies and for shorter pulses, 
which duration can be comparable with atomic natural 
lifetime. Rigorously speaking the EIT scheme can be 
also applied for storing of the short pulses if the Rabi 
frequency is large such that ft c ^> 7. But in this case 
the slope of the dispersion curve at the EIT point is al- 
ways less steep than in vicinity of the Autler-Townes res- 
onance. Thus the delay effect would manifest itself more 
effectively for the pulse which spectrum is located closer 
to the Autler-Townes resonance point. 



B. Dynamics of the probe pulse under the Raman 
process 

The general formalism developed in the previous sec- 
tion can be illustrated by the propagation of a probe 
pulse with a Gaussian intensity profile 



I Em (t) 1 2 OC L m (t) 

(Aft) 



(2tt) 1 /2 



exp 



Kf)f 



— -Aft 2 
2 



T 

t 

2 



(4.1) 



where for the sake of convenience we rescaled the field 
amplitude to unit "pulse energy" . The input field ampli- 
tude a(t) has the following Fourier component 



an 



(27T) 



1/4 



(Aft/2) 1 / 2 



exp 



QT ft 2 



Aft 2 



(4.2) 



The time shift T/2 can either be associated with the cen- 
tral point of the control pulse or, in a general situation, 
just defines the time of the probe pulse arrival. The time 
or spectral extension of the pulse profile is described by 
its spectral variance (Aft/2) 2 . 

As an illustrative example of the A-configured transi- 
tion we consider the Zeeman states in the hyperfine man- 
ifold of cesium atoms. As an initial coherent state, which 
can be prepared by optical pumping repopulation pro- 
cesses, we consider the central Zeeman state of the lower 
(ground) hyperfine sublevel \m) = \F = 3, Mq = 0), 
which is a typical situation for experiments with cesium 
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FIG. 3: (Color online) Spectral overlap of the the probe pulse 
amplitude an with the sample susceptibility x(^) — x'(^) + 
i\' (Q.) near the frequency shifted quasi-energy component 
of the Autler-Townes doublet, see figure The overlap is 
controlled by mismatching the carrier frequencies of the probe 
and coupling modes S A = A — A. The graphs are plotted for 
the following parameters: A = —507, T = 507 , Afi = 
2tt/T, £l c — I57. The spectra are shifted such that the origin 
of the plot indicates the reference point of the EIT two photon 
resonance. 



atomic clocks. The upper state |n), see figure [U we will 
associate with the relevant Zeeman state of the lower hy- 
pcrfine sublevel in the excited state of the Di-line, such 
that \n) = \F = 3,M = -1). The state \m'), coher- 
ently coupled to \m), is the relevant Zeeman state of the 
upper (ground) hyperhne sublevel, such that for cesium 
|ra') = \Fq = 4, Mo = —2). This convention in selec- 
tion of atomic transitions lets us define the scaling factor 
(|3.7[) . but all the qualitative results, discussed below, will 
be generally valid for any A-scheme. We can also point 
out that the considered configuration of atomic transi- 
tions is well approximated by A-scheme up to detunings 
A more than hundred MHz, since the nearest hypcrfinc 
sublevel in the upper state is separated by 1168Af_ffz. 

In figure [3] we show how the spectral profile cvq over- 
laps with the real and imaginary parts of the suscepti- 
bility spectral component x(fi) = x'(^) + *x"(^) when 
the atoms are dressed by a stationary monochromatic 
coupling field. The overlap of the pulse spectrum with 
the sample susceptibility can be controlled by mismatch- 
ing the carrier frequencies/detunings of the coupling and 
probe modes 8 A — A — A, see figure [5] The origin of the 
plot in figure [3] indicates the location of the EIT point 
and for the given spectra the probe pulses do not overlap 
at this point. 

Variation of the spectral location of the probe pulse 
near the Autler-Townes resonance leads to different over- 
laps with the real and imaginary parts of the sample 
susceptibility. As a consequence, the influence of ab- 
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FIG. 4: (Color online) Time dependence of the output pulses 
with the same initial spectral characteristics as in figure [3] 
which have passed the optically thick sample with no\ 2 L = 
25. The numbers in the brackets indicate the transmittance 
efficiency. The red curve corresponds to coincidence of the 
pulse carrier frequency with the Autler-Townes resonance and 
its oscillations indicate the spectral splitting of the pulse by 
this resonance. 

sorption (incoherent scattering) and delay effects on the 
pulse propagation will be also different with variation of 
SA. The negative role of absorption essentially falls off 
with deviation of SA from the exact resonance point (red 
curve in figure [3]), but the positive tendency associated 
with pulse delay can survive. This is illustrated in fig- 
ure |4] by tracking the time dependence of the outgoing 
pulses, which had the same initial spectral characteris- 
tics as in figure [3] and passed through the optically thick 
sample with n X 2 L = 25. As follows from the plotted 
dependencies, when absorption is more or less negligi- 
ble and the transmittance efficiency is around 70%, the 
pulse delay is still a perceptible effect. Let us point out 
that this predicted transmittance efficiency for delayed 
pulse is a certain indicator that the Raman protocol in 
the discussed conditions, coordinated with the currently 
existing experimental capabilities, can be applied, for ex- 
ample, for quantum storage of a single photon state with 
unknown polarization [l2j ■ In this situation it overcomes 
the best classical estimate "2/3" , see Rcf.[27] for detail. 
It is also noteworthy that in the perfect resonance situa- 
tion the outgoing pulse splits into a double-bell spectral 
profile, which in turn leads to oscillating time dependence 
of intensity, see figure [4] 

C. The pulse storage and retrieval 

In figure [5] we demonstrate the complete memory 
protocol, when the coupling light is switched off and 
switched on again after a controllable delay. The input 
pulse characteristics and the sample depth are given by 
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the same values as in figures [3] and |U From the plotted 
graphs it is evident that after a round of fast transient 
processes (poorly resolved in time scale of these graphs) 
the retrieved parts of the pulses perfectly reproduce the 
delayed parts of the original pulses. This can be seen by 
direct comparison with figure 01 The transient oscilla- 
tions have a frequency around | A| ~ 5O7 and a relaxation 
rate around 7, such that they extend on very short time 
compared with the pulse duration. For storing such a 
long pulse, in the considered configuration, the efficiency 
around 15% is achievable for effective overlap with the 
dispersion part of the sample susceptibility, see figure [3J 



Let us consider now how the Raman protocol can be 
designed for storing shorter pulses with duration compa- 
rable to the atomic natural lifetime. In figure [5] we show 
the time dependence for a set of output probe pulses, 
when the coupling field is applied exactly at resonance 
to the atomic transition, such that A = 0. In this case 
the spectral widths of both components of the Autler- 
Townes doublet are around the natural decay rate and 
the Raman scheme can be applicable for delay of short 
probe pulses. For the graphs plotted in figure [H] the ar- 
rival time is given by T/2 = 7 _1 and for Ail — 2tt/T the 
duration of the pulse, described by the time profile (|4.ip . 
is approximately the same as the atomic natural lifetime. 
Other relations and parameters are specified the same as 
for the graphs of figures [3] and 2] As a consequence the 
time dependence for the transmitted pulses, shown in 
figure El for varied mismatching of their carrier frequen- 
cies from the reference frequency of the coupling mode, 
look qualitatively similar to the relevant time dependence 
shown in figure [4] However the quantitative difference is 
great since the pulse duration and delay time are much 
shorter in the case of figure [B] The qualitative similarity 
in propagation of long and short pulses under stimulated 
Raman protocol is explained by the similarity of the ab- 
sorption/dispersion spectral profile of the Autler-Townes 
resonance for |A| ^S> £l c ^> 7 (figHJ) and for A = 0, 
«c » 7 (fig©- 



The situation changes dramatically if we look at the 
complete memory protocol applied to short pulses, see 
figure[71 As seen from this graph the role of transient pro- 
cesses becomes much more important in this case. There 
are evident manifestations of strong Rabi oscillations af- 
ter switching off the control pulse as well as after switch- 
ing it on again. These Rabi oscillations are damped on 
a time scale associated with the atomic natural lifetime. 
This time scale coincides with the pulse duration itself, 
such that the Rabi oscillations always accompany the re- 
trieval step of the memory protocol under the described 
conditions. Let us point out that the calculations pre- 
sented in figure [7] were done for a realistic value of the 
atomic optical depth and are feasible for experimental 
verification. 





SA=-0.95y (0.67/0.12) 




8A=-1.0y (0.47/0.15) 




8A=-l.ly (0.15/0.10) 




8A=-1.2y (0.4/0.12) 




5A=-1.25y (0.5S/0.098) 


/ Av 
/ // \ 

/ ' 


\ 





0.0 0.5 1.0 1.5 2.0 2.5 

t/T 



FIG. 5: (Color online) Time dependence of the probe pulse 
propagating through the sample when the coupling light is 
switched off and switched on again after a controllable delay. 
The input pulse characteristics and the sample depth are the 
same as in figures [3] and [4] The numbers in brackets indi- 
cate the transmittance/retrieval efficiency with respect to the 
input pulse. 



V. CONCLUSION 

The presented results and our discussion did not pur- 
sue the precise optimization of the pulse characteristics 
for efficient retrieval. But we can make some general ob- 
servations, concerning the applicability of the quantum 
memory protocol. Depending on the type of the memory 
protocol it may be or may not be sensitive to the fre- 
quency detuning of the coupling mode from the unper- 
turbed atomic resonance transition. For the EIT regime 
and in pure A configuration the probe light group veloc- 
ity is not sensitive to the location of the transparency 
point. If the on-resonance optical depth of the sample 
is large enough the same integral efficiency for the probe 
pulse storage and retrieval can be probably achieved at 
any detuning by proper adjustment of the original pulse 
parameters. In free space configuration the probe pulse 
stored under EIT protocol should preferably have dura- 
tion much longer than the atomic natural lifetime. 

In contrast for the Raman protocol the pulse dynam- 
ics is differently described under resonance and non- 
resonance conditions. An important consequence of our 
analysis is that the Raman protocol can be applicable 
for storing the pulses with duration comparable to the 
atomic natural lifetime. The main experimental difficulty 
in achieving the natural lifetime limit is the large mag- 
nitude of the Rabi frequency of the coupling field, which 
in such a case would be hard to stabilize. As we have 
pointed out and supported by our numerical simulations 
presented in section IIV1 the EIT and Raman protocols 
do not compete but actually complement one another. 
Each protocol has certain advantages in the area of its 
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FIG. 6: (Color online) Time dependence of the output pulses, 
when the coupling field is on resonance with the atomic tran- 
sition, A = 0. The arrival time of the probe pulse is given by 
T/2 = 7 _1 and other parameters are the same as in figures 
|3]andH An = 2-k/T fi c = 157, n X 2 L = 25. The spectral 
overlap of the probe pulse with the lower frequency compo- 
nent of the Autler-Townes doublet is controlled by mismatch- 
ing the carrier frequencies of the probe and coupling modes 
S A = A. The numbers in brackets indicate the transmit- 
tance efficiency. The red curve corresponds to coincidence of 
the pulse carrier frequency with the Autler-Townes resonance 
and its oscillations indicate the spectral splitting of the pulse 
by this resonance. 



applicability. From the practical point of view for both 
the protocols the retrieval efficiency is expected to be 
high enough for sample with on-resonance optical depth 
around hundreds and at such depths would be quite tol- 
erant to the variation of the probe pulse shape. 

We have also considered the manifestation of transient 
processes associated with switching off /on the control 
pulse. As shown by our numerical simulations the stor- 
age of short coherent pulses via the Raman protocol is 
accompanied by Rabi-type modulations of the retrieved 
pulse. The decay time of the Rabi-type oscillations and 
the pulse duration have approximately the same order of 
magnitude in this case. In the context of the memory 
protocol this is not a noteworthy effect when the coher- 
ent probe pulse approaches the single photon state and 
the encoded quantum information is associated with the 
polarization degrees of freedom of the probe light, see 
Ref. [T2| • But the effect becomes very important for con- 
tinuous variable schemes because it essentially modifies 
the spectral properties of the probe light as a carrier of 
quantum information. In this sense we point out that un- 
desirable manifestation of the transient processes is usu- 
ally less important for pulses, with duration much longer 
than the atomic natural lifetime. If storage and retrieval 
of such long pulses is controlled by the coupling field ap- 
plied in the wing of the atomic resonance line, then the 
modulation frequency of the Rabi-type oscillations is in- 



FIG. 7: (Color online) Time dependence of the probe pulse 
propagating through the sample when the coupling light was 
switched off and switched on again after controllable delay. 
The input pulse characteristics and the sample depth are the 
same as in figure [6] The numbers in the brackets indicate 
the transmittance/retrieval efficiency in respect to the input 
pulse. The retrieved parts of the the probe pulses have evident 
manifestation of the Rabi-type oscillations associated with the 
coupling field. 



creased and the degradation of the transient processes 
has much shorter time extension than the pulse duration 
itself. 

The above recommendations are quite general and are 
not constrained by unique properties of the A-scheme. 
Let us briefly comment how the developed approach can 
be further generalized to the situation of real alkali atom 
with attention payed to the hyperfine structure in the up- 
per state. It is important to take into consideration the 
hyperfine effects when probing the system in the wings of 
either the Di or D2 lines. For the coupling field tuned ex- 
actly in atomic resonance or detuned from it by a distance 
of few 7 the A-type approximation seems quite realistic 
and our results are applicable. For the practically impor- 
tant situation when the Rabi frequency is expected to be 
much less than the hyperfine splitting (Sl c <C AE\ lp i) 
and in the case of far detuning it is natural to expect 
a generalization of the Heisenberg approach introduced 
in Ref.JlTj by adding the scattering terms and the re- 
spective Langevin forces in. The important point for in- 
cluding the losses is that the incoherent scattering should 
be considered not only from the atom-associated compo- 
nents of the Autler-Townes structure (which in the limit 
fi c >7 <C |A| reproduce the atomic resonances) but also 
from the field-associated component (which in the same 
limit is located near the field mode), see figure It is 
expected that the field-associated component will be sig- 
nificantly modified due to effects of hyperfine interaction. 
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APPENDIX A: THE ATOMIC GREEN'S 
FUNCTIONS AND SAMPLE SUSCEPTIBILITY 

1. Monochromatic coupling field 

The retarded-type Green's function of the excited 
atomic state dressed by interaction with the vacuum 
modes is expressed by the standard Dyson equation 



n 

- 



. (A1) 

and then interaction with the coupling field modifies it 
to the following form 



i 

■ t 



(A2) 

The thin solid and wavy lines are respectively non- 
perturbed vacuum Green's functions of the atom and 
field. Other notation is explained in the main text of 
the paper. 

For excitation of the atom by a monochromatic cou- 
pling mode these equations can be solved by Fourier 
transformation. In the reciprocal space one has 
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Here k and u) = Wk are respectively the wave vector and 
frequency of the coupling mode and denominator is ex- 
pressed by the two resonance components of the Autler- 
Townes doublet 
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(A4) 



where u) nm i — (E n — E m >)/h is the transition frequency 
and 7 is the rate of natural decay of the upper state. 
Applying the reverse Fourier transform over the " energy" 



The reverse Fourier transform from momentum to spatial 
coordinates 
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cannot be evaluated analytically in the general case. 
However when considering ultracold atoms with temper- 
ature significantly lower than the Doppler limit the effect 
of atomic motion is rather weak and can be taken into 
account as a small correction to the approximation of 
stationary atoms. For the Raman process, when the fre- 
quency detuning between the coupling mode and atomic 
frequency is large enough, the effect of atomic motion be- 
comes even less important in the calculation of the sam- 
ple susceptibility. But it gives the principle limitation for 
the storage time of the memory protocol as discussed in 
subsection IIII Bl and commented below in sec. I A 41 



2. Coupling with a rectangular time envelope 
profile 



If the action of the coupling mode is limited in time 
and described by a rectangular envelope profile (|3.9[) in- 
side the interval < t < T, then the Green's function 
GwnT^ (p; t, t 1 ) defined in the mixed representation for dif- 
ferent relations between times t and t' is given by 

1.1) t' < and t < 
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1.3) t' < and T < t 
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-B„ + (p,a;) - £ n _(p,w) 

£ n -(p, w) - £„+(p, w) 

2.1) < i' < T and < i < T 
G?( l --)(p;*,0 = -^(i-0 
[E n+ (p,u) - B„_(p,w) 

. ^n- (P, - E m > (p, U) 



(A9) 



■ e b 



fE„_(p,w)(t-t') 



£?„_(p,u;) - £ n+ (p,w) 
2.2) < i' < T and T < f 



(A10) 



£„+(p,Cj) - E m r(p,U}) 



£ n+ (p,u/) - E n _(p,w) 

£ n _(p,^) - £ m /(p,o;) <B„_(p lW )(T-t') 



(All) 



£ n _(p,w) - £ n+ (p,w) 
3) T < t' and T < t 
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These expressions for the Green's function in different 
time domains can be built up by clipping the partial so- 
lutions for different excitation regimes. These include 
the atomic dynamics driven by either interaction only 
with the vacuum modes (|A7[) . (|A12[) or jointly with the 
vacuum and coupling modes (|A5|) . The required bound- 
ary conditions are dictated by continuity at the points 
t = 0,T. 



3. Susceptibility of the sample in response to the 
probe field 

The susceptibility of the atomic sample in response 
to the probe field is given by the polarization operator 
contributing to the Dyson equation (|3.ip 
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where no(rj_, z) is the spatial distribution of atomic den- 
sity. The total energy of the atom in the ground state 

E m (p±,p z -hk)=E m + ^+ (Pz ~ %l)2 (A14) 

2m 2m 



is described by the kinetic part, written as function of 
atomic momentum p = (pj_,p z ) of the excited state and 
shifted by the momentum of the photon from the probe 
mode hk. Such a momentum conservation in the coherent 
process (|3.ip is also taken into account in the argument 
of the Maxwell distribution fo(p±>Pz — hk). 

It is convenient to define the slow-varying component 
of the susceptibility 



(A15) 



which describes the dynamics of the slow- varying ampli- 
tude of the probe field, see Eqs.O, l[3To]) . <\5T2\i , (|3~T3)) . 



4. The Green's function and sample susceptibility 
modified by interruption of the control pulse 

In the general case the retarded-type Green's function 
responsible for atomic dynamics, dressed by two pulses 
of the coupling field (see diagram p. lip ), has a rather 
cumbersome analytical structure. The result can be sim- 
plified for the most important situation when time t' be- 
longs to the first control pulse and time t belongs to the 
second pulse. Let us denote the delay between the two 
pulses by Tjj and superscribe all parameters associated 
with the second pulse, such as duration, Rabi frequency, 
transition matrix elements etc., by a prime. Then for 
< t' < T and T + T D < t < T + T D + T' the Green's 
function, expressed by diagram (|3.1ip . is given by 



G { n jJ(p;t,t') 



V ,v* , 

7im' nm' 



{ E 'n+ - E 'n-)( E n+ - E n -) 
, e Ti E ' n+ (t-T-T D ) _ e -^E' n ^(t-T-TD)] e -iE m ,(p,u,)T D 



-|E„+(T-t') _ ±E„-(T-t') 



(A16) 



were the matrtix element V^ m , is related to the retrieval 
coupling field. The susceptibility component responsible 
for the recovering process, described by the first diagram 
on the right hand side of Eq. (|3.10p , is given by 

1 



X//(r±,z;t,t') = - 

d 3 p 
(2irhy 



de -)nml 

no(r±,z)f (p±,p z - hk) 



xe fi 



Up^-fiS)(t'-t) G e-) (p . t) ^ (A17) 



and its slow-varying component x//(t±, z\ ^ 0> con- 
tributing to Eq.(3.12), is defined by the transformation 

(|ATo)) . 

It is important to recognize that validity of expression 
(|A17|) is restricted by the global effect of atomic motion 
for the entire memory protocol, which includes both the 
storage and retrieval stages. In our derivation of the sus- 
ceptibility (|A17|) . driven by two control pulses, we have 
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assumed that during the delay time Td atoms do not tial inhomogeneity associated with their original density 
noticeably change their location on the scale of the spa- distribution n (r±, z). 
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